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Abstract. We report on the theoretical description of secondary electron emission due to resonant charge 
transfer occurring during the collision of metastable N2( 3 X'^") molecules with dielectric surfaces. The 
emission is described as a two step process consisting of electron capture to form an intermediate shape 
resonance N^~( 2 77 g ) and subsequent electron emission by decay of this ion, either due to its natural life time 
or its interaction with the surface. The electron capture is modeled using the Keldysh Green's function 
technique and the negative ion decay is described by a combination of the Keldysh technique and a rate 
equation approach. We find the resonant capture of electrons to be very efficient and the natural decay to be 
clearly dominating over the surface-induced decay. Secondary electron emission coefficients are calculated 
for AI2O3, MgO, Si02, and diamond at several kinetic energies of the projectile. With the exception of 
MgO the coefficients turn out to be of the order of 10 _1 over the whole range of kinetic energies. This rather 
large value is a direct consequence of the shape resonance acting as a relay state for electron emission. 

PACS. 34.35.+a Gas-surface interactions - 34.70.+e Charge transfer - 79.20.Hx Secondary electron 
emission 



1 Introduction 

Secondary electron emission from dielectric surfaces due 
to various projectiles is a fundamental process of great 
technological significance. The associated electron yield is 
usually represented by the secondary electron emission co- 
efficient 7 e which denotes the average number of electrons 
released by a single collision of a particular type. Thus, j e 
is specific to the surface material, the projectile, and the 
particular collision process under consideration. 

To illustrate the importance of secondary electron emis- 
sion from dielectric materials we first consider an atmo- 
spheric pressure dielectric barrier discharge (DBD). Var- 
ious technological processes like polymer surface treat- 
ment, thin film coating, sterilization, gas flow control and 
ozone production |1I2| are based on this type of discharge. 
All of these are sensitive to the discharge's stability and 
homogeneity. Secondary electron emission from the elec- 
trodes increases the DBD's self-sustain and is thus crucial 
to the discharge's operation mode. In particular, there is 
some evidence that the stability of DBDs is controlled by 
secondary electron emission from plasma boundaries by 
impacting metastable species Although DBDs have 
been studied for more than 20 years, their basic mecha- 
nism is still not entirely understood pQ. To overcome this 
deficiency numerical investigations are inevitable |2l3j . These 
studies require j e as an input parameter, which either has 
to be measured experimentally or calculated from first 
principles. 



Another application are plasma display panels (PDPs) , 
which are used in large size TV monitors. PDPs are usu- 
ally based on AC driven high pressure microdischarges [I] 
in noble gas mixtures [5] and often employ electrodes coat- 
ed with AI2O3 or MgO jBj. The panel's firing voltage, that 
is, the voltage required for discharging, is directly influ- 
enced by secondary electron emission from the dielectric 
electrodes [5J, for instance, due to impacting metastable 
species. In order to improve the performance of PDPs, 
theoretical investigations of the employed materials, the 
overall discharge geometry, and the operation mode are 
necessary [4] . These studies again require the emission co- 
efficient 7 e as an input parameter. 

The demand for reliable j e data for use in numerical 
simulations of PDPs, DBDs, and related gas discharges 
contrasts with the very sparse and often uncertain avail- 
able reference data. In many cases rules of the thumb 
have to be used [7] . Especially sparse is the data base for 
secondary electron emission due to metastable molecules. 
The objective of the present work is therefore to contribute 
to the filling of this particular gap. We aim to establish 
a model that allows for easy calculation of the secondary 
electron emission coefficient due to impacting metastable 
molecules on dielectric surfaces, taking N 2 ( 3 X'+) as an ex- 
ample. Other projectile-surface combinations can be mod- 
elled in the same spirit. 

The model we propose is an effective one, stripped of 
some of the microscopic details of the system and charac- 
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Fig. 1. Energetic structure of the different dielectric mate- 
rials considered in this work characterized by the electron 
affinity s a , the size of the energy gap e g , and the width of 
the valence band Aev- For SiC>2 the energetically blocked 
upper valence band is marked in shaded grey. Also in- 
dicated is the variation of the molecular hole level with 
distance. For an explanation of the molecular image shift 
see Sec. |2.1| Note that the figure is to scale in terms of 
energy units. 



terized by a small number of material specific parameters 
instead. In a previous work [5] we presented such a model 
for the case of Auger de-excitation of N2( 3 I7+) on alu- 
minum, corresponding to the reaction 



N 2 ( 3 i;+) + e fc 



N 2 (^ 



(1) 



with efc and e q denoting an electron within the surface 
and a free electron, respectively. Another important pro- 
cess capable of generating secondary electrons by impact 
of N 2 ( 3 I?+) is the combination of electron capture by res- 
onant charge transfer (RCT) to form the N^~( 2 77 s ) shape 
resonance and subsequent decay of the intermediate ion, 
corresponding to 



N 2 ( 3 r+) + e fe 



RCT 



N 2 -( 2 7I 5 



RCT 



N 2 (^ 



(2) 



Except for diamond, the Auger channel ([T]) is energetically 
suppressed for all of the dielectric materials considered in 
this work whereas the RCT channel ([2| is energetically al- 
ways possible (see Fig.[T]). Hence, we focus on the resonant 



reaction ^ which was also found to be very efficient when 
N^ 3 !^) collides with a metallic surface |9ll0j . In order 
to describe this process, we will use a generalization of the 
model presented in [8] and apply it to dielectric surfaces. 
Since the adaption is straightforward we only give a brief 
sketch of the model and refer the reader to [5] for details. 

The solid surface is treated as a potential well with a 
depth 



Vn = 



-e a — e g — Aey 



filled up to the top of the valence band 



£v 



(3) 



(4) 



Here e a , e g and Aev denote the solid's electron affinity, its 
energy gap and the width of its valence band, respectively. 
The wave functions of such a potential well are easily cal- 
culated and may be found in [8] . The particular electronic 
parameters of the dielectric materials investigated in this 
work [111121131141151161171 18 19 20 21] are summarized in 
Table [U 

The molecule is characterized by a two level system 
consisting of an initially empty ground state level and 
an initially filled excited level 1 . The molecular wave func- 
tions associated with these two levels are calculated from 
a LCAO combination of atomic nitrogen wave functions, 
identifying and 1 with the 27r u and 27r g molecular or- 
bitals [8 . The ionization energies of the two lev- 
els [22123] are listed in Table [2] for the unperturbed mol- 
ecule. As will be shown in Sec. |2.1| both of these levels 
shift downwards upon approaching the surface, provided 
that only the RCT channel ^ is considered. 

We will describe the resonant reaction pi as a two- 
step process. First, we will employ the Kelaysh Green's 
function technique to model the resonant electron capture 
into the initially vacant ground state level 0. This step of 
the process is depicted in Fig. [T] The second step consists 
of resonant electron emission by decay of the intermediate 
negative ion. The decay can be due to the ion-surface cou- 
pling, which is mediated by the image interaction, or due 
to the ions natural life time, that is, due to auto-decay. 
Considering a N^( 2 77 g ) ion colliding with the surface we 
will model the surface-induced decay again by a Keldysh 
approach, the auto-decay, however, by a rate equation ap- 
proach. Electron capture and release we will be coupled 
by a rate equation scheme that is based on a conditional 
decay reaction. 

The outline of the remainder of this article is as follows. 
In Sec. [2] we will setup the quantum mechanical model of 
the molecular electron-capture and conduct a quantum- 
kinetic calculation on top of it. The treatment of surface- 
induced decay and auto-decay as well as their coupling to 
the initial electron capture will be presented in Section [3] 
Finally, we show and discuss numerical results for selected 
surfaces in Sec. [4] and give a conclusion in Sec. [5] 



Johannes Marbach et al.: Resonant charge transfer at dielectric surfaces 





e a [eV] 


£a [eV] 


4ev [eV] 


4 


M203 


1.0 01] 


8.7 [n] 


11.8 [12] 


9.0 [13] 


diamond 


0.52 [H] 


5.48 [H] 


21.0 [T5] 


5.66 [H] 


MgO 


0.925 [T7] 


7.6 [E] 


8.5 [E] 


9.83 [19] 


SiOa 


0.9 [2D] 


13.8 [21] 


5.3 [5T] 


3.9 hh] 



Table 1. Electronic parameters for the dielectric mate- 
rials considered in this work: electron affinity e a , energy 

irn Innnn Kfinrl /li-. . nnrl o+- r\ 4- n r* -vr\ 1 ri 
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gap £ g , valence band width Aey, and static relative di- 



electric constant e\. The electron affinity of diamond and 
MgO is calculated from a mean value due to lack of con- 
sistent data. Note, for SiC>2 only the lower valence band is 
considered. The upper valence band starts at 8.9 eV and 
has a width of 2.7 eV [21]. It is thus energetically well 
separated from the molecular vacancy level £0 and does 
therefore not contribute to the electron capture. 



2 Resonant electron capture 
2.1 Model 

The molecule-surface system is described by the following 
Hamiltonian 

h 

c fc c c k c o ) > 



(5) 



+ 



where k denotes electronic states within the solid's valence 
band and labels the molecule's ground state level, which 
is initially empty. Since we are considering the capture of 
an additional electron by the metastable state N2( 3 i7+), 
the level equals the lower ionization level of N^"( 2 77 g ). 

To understand the variation of £ with time we fol- 
low [24J and consider a negative nitrogen ion situated at 
a distance zr from a dielectric surface, the relative static 
permittivity of which we denote by e r . Since the ion is neg- 
atively charged its overall energy is lowered by the value 
of the image potential at zr, for which we use the classical 
expression (see however below) 



Vi{z R ) 



4-1 e 2 1 

e b r + 1 47T£ Az R 



(6) 



In order to obtain the effective ionization energy eq(zr), 
we virtually move the negative ion from the position zr 
to a point infinitely far away from the surface, ionize the 
electron bound in the state to the vacuum level, and re- 
locate the resulting neutral molecule back to the original 
distance zr. The net energy change during this process 
determines the effective ionization energy 



e (z R ) = -Ae = -[ -V^zr) + (-eg°) 
= eo+Vi(z R ) . 



N 2 -yz R 



(7a) 
(7b) 
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Table 2. Ionization energies of the molecular ground 
state level and the excited level 1 for a free nitrogen 
molecule in its ground state N 2 ( 1 i7+) and its metastable 
state N2( 3 i7+) as well as for the negative ion shape reso- 
nance N^( 2 J7 g ). 



Since this procedure is independent of the specific molec- 
ular level, we conclude that upon approaching the surface 
each and any of the negative ion's ionization levels shifts 
downwards. The reader may note that the decrease of the 
ionization energies does not correspond to an actual en- 
ergy change of the negative ion. The ionization levels are 
defined solely in terms of the removal of an electron from 
the molecule and, thus, do not imply an energy change for 
a bound electron. 

Employing the usual trajectory approximation [25] to- 
gether with the trajectory 



R(t) = z R (t) e z = (v \t\ + z ) e 2 



(8) 



the relation (7b I can also be understood as £q(£). 

Generalizing Gadzuk's j26| treatment of mono-nuclear 
projectiles to bi-nuclear projectiles, the tunnelling matrix 
element is given by 

V u (t)= Jdrn{r)[V M (r,t) + V S M(z)]#o(r v (t)), (9) 

where the time dependence arises from the molecule's mo- 
tion and r v equals the vector r as seen from the molecular 
reference frame, that is, 



rJt)=f2(<p)[r-R(t)} , 



(10) 



with an appropriate rotation matrix f2(ip) [S]. 

The expression Vm specifies the Coulomb potential of 
the screened molecular nuclei 126 



V M (r,t) 



Z cS e 2 
47T£o e r {z) 



1 



1 

H 2 (i)| 



with Z e g ss 4 [8] and 



R 1/2 (t) = R± 



(11) 



(12) 



denoting the positions of the two nuclei, which are given 
in terms of the center of mass coordinate R and the inter- 
nuclear vector q. The static dielectric constant e r {z) is 
approximated by 



e r (z) 




(13) 
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The term Vsm labels the electron's self image interac- 
tion potential 26J 

Vsm(z) = 0(z) [G(z c ~ z)V + 9(z ~ z c )Vi(z)] . (14) 

Here we have used the usual truncation of the classical 
image potential close to the surface [57] with 



and ctq given by the infinite series 



<r (t,t') = 1- / dfi / dt 2 A (h,t 2 ) 



dti / dt 



dt? 



* 3 (22) 

dU A)(*l,*2)A)(t3,«4) 



Vi(z c ) = V 



(15) 



The wave functions of the metal electron Wk and the 
one of the molecular electron <Fq are taken from [8]. We 
treat the magnetic quantum number to, which is intro- 
duced within the LCAO representation of the molecular 
wave functions [8], as an external parameter. 



2.2 Quantum kinetics 

The Hamiltonian ^ describes a standard problem of many 
particle physics. It can be treated exactly using the Keldysh 
technique [55] . The self-energy solely consists of a first or- 
der retarded and advanced term 



The quantity Aq emerges from the self-energy and reads 
A (h,t 2 ) = pXV fc *(fi)V fc (i2) 



k (23) 



By employing now the Dyson equation for Gq\, with 
the self energy ( 16b), the occupancy of the molecular level 



( 18b ) can be written as 

rt 



n (t) 



dii / dt 2 A>(ti,*2)<ro(f,*iK(M2) • (24) 



yR/A, 
^kO l*!'^ 



yR/A, 

^ok {piih 



(16a) 
(16b) 



Here Z\o is equal to Aq with an additional factor nk{to) un- 
der the fc sum. In gen eral this subtlety would deny further 
simplification of Eq. (24 1. However, since for a dielectric 



surface the active energy band is the valence band, which 
is completely filled, the terms Aq and Aq are essentially 
The Keldysh component of the self-energy vanishes, the same. Note that even for metallic surfaces with an 
the Dyson equation for the Keldysh Green's function can incompletely filled conduction band these terms will still 

be almost equal, as long as the molecular resonance stays 



thus be solved to give 

G^(M') = -i£[l-2n 7 (t )] 

7 



inside the occupied 



portion of the band. 



(17) 



To simplify Eq. (24), we first note that when t is kept 



constant the derivative of Eq. ( 22 ) with respect to t' can 
be shown to read 



where a and /? can take the values or k and 7 runs 
over all states {0, fc} of the system. Note that to is short- 
hand for —00. Equation (17) can be used to calculate the 
occupation of the molecular level 

1 



d<T (t,t) 
dt' 



df<7o(t,f)A>(t,0 



Using this expression we can now rewrite ( 24 ) into 



no(t) = -[l-iG&(t,t)] 

= 5> fe (t )|G R fe (Mo)| i 

k 

where we have used the relation [28129] 

/J Gg-y (*) *o) Gya (frh f ) = 1 



(18a) 
(18b) 

(19) 



n (t) = 25R 



dt 2 cro(t,ts 



d<7 a (t 1 t 2 ) 
dU 



2[l-o- (t,toK(t,t )] 

A<y* {tM) 



23? 



dt 2 cr (t,t 2 )- 



ta 



dU 



(25) 

(26a) 
(26b) 



where we have used integration by parts in the last line. 
Combining (26a) and (26b), we obtain 



Evaluation of the Green's function G^ fe appearing in ( 18b) 
requires knowledge of the diagonal component Gqq, the 
Dyson equation of which can be solved iteratively to yield 

G*(t,t') = G*; \t,t')a (t,t'), (20) 

with 



n (t) = 1 - oo(Mo)<>o(Mo) 



(27) 



We have not yet used any approximation to ctq- The 
result (27) is thus exact. Since Vk is small it is however 
sufficient to approximately re-sum the infinite series ( 22 ) . 
Using the exponential re-summation technique described 
in [8] and retaining only the lowest order term we obtain 



R(o\ 

00 



(i,f') = -i0{t-t')e-*ti dtl eo(tl) , 



(21) 



(28) 
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which after insertion into (27) yields 

no(t) « 1 - e _ J"*o dtl /*o d * 2 ^(ti.fa) 



(29) 



Evaluating (29) at t = i we see that the initial condi- 



tion no (to) = is fulfilled. Taking the limit L — ¥ oo within 
the metallic wave function [8] , the k sum in ( 23 ) turns into 
an integral. For a numerical treatment it is convenient to 
exchange the k and Uj2 integrations in ( 29 1 . The integral 



in the exponent of ( 29 1 then takes the form 



(30) 



with 



A (k,t) = f dii y fc (t 1 )e* £fctl es-''o ldt2£o(t2) , (31) 



where the L 3 factor in ( 30 1 cancels with the L 3 / 2 factor 
in V k . 

2.3 Matrix element separation 

In order to make a numerical calculation of the ground 
state level's occupancy feasible, we need to introduce ap- 
proximations to the matrix element. 

Numerical evaluation of the full matrix element (|9| 
suggests that Vm is clearly dominating over Vsm and that 
in addition the matrix element may be neglected in the 
half space z < 0. Making use of these findings we can sep- 
arate the matrix element's time dependence by performing 
the coordinate transformation 



(32) 



which is the inverse of (10). Utilizing the relation 



\r-R 1/a (t)\ = \T-R(t)T%&(tp), 
= \&{v)KT~e*]\ 



and approximating 0(zL + z R (t)) by 0(z' v 
we then obtain 

where T kz is defined in [jj], 
mN k Z cS e 2 



, and 

V k ldr v G(z'-z )e-"<°* z '* 



x e 



e '*> +e ' 



(x<p + imy v ) 



] (33a) 
(33b) 
(33c) 

zq) as in [8] 

(34) 
(35a) 



(35b) 



Here we have rescaled r v , k and t according to 



1 



k — » xk 



1 

t -> 1 

XV 



(36) 



with x = . Furthermore, the expressions are short- 



hand notation for 



(37) 



Inserting (34) into Eq. (31) yields 



A(k,t) 



dt\ e 



x e n^ /o ldt 2 eo(«2) 



(38) 



which enables us to calculate the involved space and time 
integrals independently within a numerical simulation. 



3 Decay of the negative ion 
3.1 Surface-induced decay 

The surface-induced decay reaction can be described by 
the Hamiltonian 



(v q (t)cl Cl + V*(t)c q c\) 



(39) 



which is structurally identical to the resonant capture 
Hamiltonian ([5| with the associations k «-» q and Ool. 
In (39) E\ is the upper ionization level of the negative 



nitrogen ion and e q denotes the energy of the free elec- 
tron states, which are centered around the moving mol- 
ecule. Both of these energy levels shift downwards upon 
approaching the surface, as can be shown using the anal- 
ysis presented in Sec. |2.1| 



(40a) 
(40b) 



e 1 (t)=sT + V i (z R (t)), 
e q {t)=e q ° +Vi{z R {t)) . 



This fact implies that the difference £i(t) — £ q (i) is con- 
stant in time. 

The tunnelling matrix element V q is given by 



V q (t) = a(e™,t) dr^ q (r v (t))^(r v (t)) 



x [Vs(z) + V SM (z)] , 
with tZ'i as defined in [H] and & q given by 

V q {r) = Xe*** . 

L2 



(41) 



(42) 



Since the electron is emitted from the moving mole- 
cule, we have to evaluate in Eq. (41 ) the wave function of 
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the free electron, W q (r) : in terms of the molecular coordi- 
nates r v and the wave vectors q v which are the electron's 
wave vectors as seen from the molecule's reference frame. 



The plane wave ansatz (42) is valid in the present case, 



because the emitted electron leaves a neutral molecule be- 



hind. Conceptually, using W qxp [r v ) in Eq. (41 ) corresponds 
therefore to the two-center Coulomb wave used in [5] but 
with zero effective nucleus charge. 

In (41 1 we introduced the term cv(e^,i) which ac- 
counts for the trapping of slow electrons in the image 
potential close to the surface. When calculating the occu- 
pancies of the molecular levels, a must equal unity since 
the decay of the molecular state does not depend on the 
emitted electron being able to escape. On the other hand, 
when calculating the number of emitted electrons, a needs 
to cut off the matrix element when the overall energy of 
the emitted electron is smaller or equal to zero. 

The overall energy of the electron is given by the sum of 
the electron's perpendicular energy and the image po- 
tential at the position of emission, which may be approxi- 
mated by the molecule's center of mass coordinate [8j . We 
thus define 



1 , mol. level 

6>(e~ + Vi(z R (t))) , q states 



(43) 



The interaction potential Vs is given through 

V s (z) = 0(-z)V o , (44) 



and Vsm is defined in Eq. ( |14| ). Employing the transfor- 
mation (32) the matrix element takes the form 



V,(t) = a(e~ i) /drviP* (fv)!Pb(»v) 



x [V s {z' v + z R (t)) + V SM {z' v + z R (t))] 



(45) 



Focusing on the two principal molecule orientations 
and using the explicit form of the wave functions together 
with the calculus of Fourier transform Eq. ( 45 ) can be 
reduced to a single integral 



v q (t) 



Li 



/OO 
-oo 



(46) 



x [V s {z' v + z R (t)) + V SM {z' v + z R (t))] 
In the perpendicular geometry (P v reads 

Q 



<P±(z' ± ,Q ± ) = -(2i)V27r 



im<P\ 



(i + Q 2 )t 

\z+\* K,(\z+\y/l + Q*) 



(47) 



IzJ 3 Ks 



with K denoting the modified Bessel function of the sec- 
ond kind and 

z ± =^±f. (48) 



For the parallel orientation <P V takes the form 
<2>ll(z[|,Q||) = AiV^sin (g x |) 

>n ° C ° B f". Uil^ S (l41V^+Q' 



(1 + Q 2 P 
\z'\i 



(49) 



-Kg 



<ilvT 



"(l + Q 2 )t 

The lateral wave vectors Q± and are given by 

Q± = (q x ,q y ) = (Qcos(# x ),Qsin(^ x )) , 
Q\\ = (q y ,q x ) = (Qcos(<£||),Qsin(<2>||)) . 



(50a) 
(50b) 



Since the Hamiltonian (|39j) is structurally identical 
to ([5]), the quantum-kinetic calculation outlined in Sec. 2.2 
holds for the current case as well. The only difference is 
that the ion's motion does not start at t = to = —00, but 
at some later time t g , at which the ion is generated by 
electron capture. This can, however, easily be accounted 
for by replacing to with t g in the previous calculation. 

Using an adapted version of Eq. (18b) we thus find for 



the occupancy of a free electron state with wave vector q 



1 



(t) = -[l-iG^ q (t,t)]=\G^(t,t s 



(51) 



Employing the Dyson equation for G ql (t,t g ) and the 
results of Sec. I2.2l the total number of emitted electrons is 
then given by 



(52a) 



f dti f dt 2 A 1 {t 1 ,t 2 )al(t 1 ,t g )a 1 (t 2 ,t g ) , (52b) 
Jt„ Jt„ 



1 - <7i(t,t g )al(t,t g ) 



with 



4i(*i>*2) = ^X>q(*i)W) 



x e 



<3 



(52c) 



(53) 



-Ef)(tl-i2) 



and ex 1 defined by an expressions analogous to Eq. ( 22 ) 
but with Ai instead of A . 



The equality of ( 52b ) and ( 52c I , although not obvious 



can easily be confirmed. To do so, we denote these two 
equations by Ci(^) an d respectively, and calculate 

their derivatives with respect to t 



dCi(*) 
d< 



dt 2 A 1 {t,t 2 )'yl{t,t g )a l {t 2 ,t g ) 



+ / dti A^t^t) alit^tg) a^t.tg) 



(54a) 



'^'■' ] (--.410 



df 



dCaft) 
At 



(54c) 
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In addition, by inspection of Eq. ( 52 ) we see that Ci (*) 
and Ca(*) both vanish at t = t g . Hence, (i(t) and CaC*) are 
solutions of one and the same ordinary first order differen- 
tial equation with one and the same initial condition. By 
virtue of the uniqueness of such a solution Ci(i) and C2W 
are thus identical. 

The time integration of the modulated matrix element 
involved in ( 52c I , requires evaluation of 



1 



1 f°° ■ 1 

VL2 J-OO 



(55) 



where we have used the variable scaling ( 36 1 and intro- 
duced the quantity 

I(t,t g ,z v )= [ dt x a(e~ tOes^S - 6 ")** 



Jt„ (56) 
x [V s (z' v + z R (t 1 )) + Vsm(^ + zn(h))] • 



Although tedious, the time integration within Eq. (56 1 
can be done analytically. For shortness, we merely give 
the result 



3-(t> tg , Zip) — 



2Vc 



x sm a 



a 

h-ti 
2 



0{t 2 - h) 0(-t g ) e* 



, «l+«2 



ia f 3+*4 . 

x e 2 sin I a 



+ 9(t)0(t 4 -t 3 ) 
t 4 -t 3 



167T£o 



0(t 5 -t g )0(-t g )e™ 

x {Ei(-ia(«-t B )) + J(H-*sgn(a))} 
- 0{t) 0(t - t 6 ) e~ iai 

x |Ei(ia(z + f)) - Ei(ia(z + f 6 )) } 
where we have introduced the following abbreviations 



(57) 





= max(i 9 ,z^ + z - Z c ) , 


(58a) 


h 


= min(0,^,t) , 


(58b) 


*3 


= max(i ff ,£ 9z ,0) , 


(58c) 


*4 


= min(z c - z' 9 - Z ,t) , 


(58d) 


*5 


/ 

— z c. — — Z , 


(58e) 


is 


= min(-i 5 ,-t gs ,,0,t) , 


(58f) 




= max(i ff) 0,i 5 ,i g J , 


(58g) 




4 - 1 ne 2 1 


(58h) 




£^ + 1 167T£ £q z 




C „ Ci 





Equations ( 57 ) and ( 58 ) hold when calculating the num- 
ber of emitted electrons n(t). They can, however, also be 
applied to calculate the occupancy of the molecular levels 
if one sets t qz equal to 0. 

In applying the analysis of Sec. |2.2| to the excited 
molecular level, we find its occupancy to be 



ni(i) = 1 - n(t) = <T 1 (t,tg)al(t,t g ) 



(59) 



In order to obtain the rate of the surface-induced decay 
process /g, we use the lowest order approximation for a± 
and take the derivative of Eq. ( 59 ) with respect to t, 



dm(t) 

dt 



23? 



<rj(t,i fi 



dai(t,t g ) 



dt 



dt 2 s R[A 1 (t,t 2 )}a 1 (t,tg)a* 1 (t,t g ) 



r s {t,t a )m{t) 



with 



r s (t,t g 



2 / di 2 R[^i(*,t 2 



(60) 



Note that in addition to the time of observation t, the de- 
cay rate r$ also depends on the time of generation t g . The 
coupling of the surface-induced decay process with the 



prior negative ion generation will be outlined in Sec. 3.3 
Before, however, we turn our attention to the natural de- 
cay of the shape resonance N^~( 2 77 g ). 

3.2 Natural decay 

The N^( 2 77 ff ) shape resonance has a mean life time due 
to natural decay of r = 1.6 • 10~ 15 s |30J, which corre- 
sponds to a decay rate of Tjv = r _1 = 0.625 • 10 15 /s. In 
earlier work [9] the line shape of auto-decay processes has 
been modelled by means of a Lorentzian distribution, also 
known as Breit-Wigner distribution. In accordance with 
this approach we define therefore the natural decay rate's 
energy spectrum f?jv(£^°,i) by 



r N (t)= / dE~ ejv (e~ f) 
Jo 

As r°° n [ 2lx 



(61) 



q(gg° cos 2 (t9 g ),t) 

x K-s^ + ^f ■ 

Here As is the line's full width at half maximum, which 
can be estimated using Heisenberg's uncertainty relation 

h 



As > 



(62) 



In (61) we again introduced the image potential trap- 



ping factor a(sq z ,t). Using its explicit form (431 we can 
carry out the angle integration in (61) and obtain 
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with 



mol. level 



' J«te®l), q states ■ (64) 



The integration boundary e a appearing in ( 63 1 is given by 

(65) 



, mol. level 

£i(zn(t)) , q states 



The rate equations for the complete electron emission 
process needs to take into account not only the decay pro- 
cess, be it natural or surface-induced, but also the preced- 
ing generation of the negative ion. It is described next. 



3.3 Consecutive decay approach 

Because of the lack of an effective Hamiltonian for the ac- 
tive electrons describing simultaneously the formation of 
the negative ion due to electron capture and the decay 
of the negative ion with subsequent electron emission we 
cannot employ the Keldysh technique to calculate the sec- 
ondary electron emission coefficient. As an alternative, we 
use rate equations with rates obtained from the Keldysh 
approach. 

The reaction chain consists of two consecutive pro- 
cesses. Considering only the molecular states we can write 



N 2 ( 3 ^+ 



A(t) 



> N 2 ( x i7+) , (66) 



where A(0 and A(i) are the reaction rates of electron 
capture and release, respectively. The scheme ( 66 1 is easily 



translated into a system of ordinary differential equations 



dn*(t) 

dt 
dn_(i) 

dt 
dn g (t) 

dt 



= -r (f)n,(t) , (67a) 
= -A(t)n_(t) + r (t)n,(t), (67b) 



= A(t)n_(f) 



(67c) 



where n*, n_ and n g denote the respective fractions of 
metastable molecules, negative ions, and ground state mol- 
ecules, which obey the relation 



n*(t) + n-(t) + n g (t) = 1 V< 
and satisfy the initial conditions 

n,(to) = l, n_(to) = 0, n g (t ) = 0. (69) 
The system ( [67] ) can be solved straightforwardly to give 

n»(i) = e- f to dtl , (70a) 

«_(*) = - fdt, dn * ih) Jh ^ (70b) 

J to d *l 



n g (t) = 1 -n*(t) - n_(f) 



The generation of negative ions is realized by filling 
the hole in the ground state level 0, that is, the lower ion- 
ization level of ( 2 n g ). The ion then decays by emitting 
an electron from the excited level 1 which is the higher 
ionization level of N^"( 2 ./7 ff ). Hence, the molecular frac- 
tions n^/^/g can be related to the single level occupan- 
cies n /! by means of 



n,(t) = l-no(t) 
n g (t) = 1 - m(t) 



n-(t) = n (t) + n 1 (t)-l . (71) 



To calculate the total number of emitted electrons n(t) 
and their energy spectrum we first note that due to par- 
ticle conservation n(t) must equal the fraction of ground 
state molecules n g (t). Hence, we can employ Eq. (67c) 



to calculate n(t). In analogy to (61) we write the decay 
rate A(*) as 



A(*) 



(72) 



where we introduced the rate's energy spectrum gi(e^°, t) 
which needs to account for trapping of low energy elec- 
trons in front of the surface due to the image potential. It 
will be given below. 



After employing ( 72 ) in ( 67c ) we directly integrate the 



resulting equation and obtain 



n(t) 



de° 



(73) 



Inserting ri_ from Eq. ( 70b ) we can use ( 73 ) to calculate 



the overall number of emitted elec trons as well as their 
energy spectrum. For t — oo Eq. ( 73 ) defines the secondary 



electron emission coefficient 7 e = n(oo). 

Finally, we need to build a link between the rates A/i 
and the previously calculated quantities A/JV- The rate 
of electron capture A is defined by 



d(l-no(*)) 



dt 



-r (t) (1 - n (t)) 



(74) 



and can be calculated from (29). The rate of negative ion 



decay A contains contributions from the surface-induced 
decay and the natural decay. For the latter the rate Av 
was explicitly given in Sec. |3.2| For the former we need to 
keep in mind that the rate rs(t,t g ) given in Sec. 3.1 also 



depends on the time of negative ion creation t g . To obtain 
the overall rate A?(*)j we therefore need to weight the 
two-time rate rs(t,t g ) with the negative ion's production 
term along the trajectory, that is, 



A(i) 



dt g r s (t,t g ) r {t g )n*(t g ) 



(75) 



to 



Employing the explicit form of rg(t,t g ) given by Eq. (60) 
together with partial integration we obtain 



(70c) r s (t)= I dti2R[4 1 (t,t 1 )](l-n,(ii)) 

Jt 



(76) 
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Fig. 2. Time evolution of the occupation of the molecu- 
lar ground state level no in front of an AI2O3 surface for 
parallel (solid line) and perpendicular (dashed line) mol- 
ecule orientation. The molecule's kinetic energy was fixed 
to 50 meV. The curves were calculated from Eq. (29). 



The total rate of decay is given by 

r 1 (t) = r s (t) + r N (t) 



(77) 



and the total rate's energy distribution gi(e^,t), which 



was introduced in (72 1, is related to the previous calcula- 
tion by means of ( 63 ) and (|76| , respectively 



To be specific we note that for constructing Q\ (e^° , t) 
the q sum involved in A\ (see Eq. ( 53 )) has to be converted 



into an integral by letting L — > 00. Afterwards one can 
integrate the angle coordinates of the q vector to obtain 
the energy spectrum of the rate for the surface-induced 
decay. The final expression for Qi (e^° , t) reads 



Q 1 (e™,t) = Q S (s™,t) + Q N (s™,t) 



(78) 



with 



0S = 



2m P 



l 3 r* /" f [ 2 * 

dti / dV q / d(p q 



h 2 K 2 J h 2 
x k jVq(*) v*{ti) e w( £ ? > - £ r > )( t - t i)| 

x ^/i" sin(i?,)(l-n«(ii)) 



(79) 



and Qx(e™,t) implicitly denned in (61) 



4 Results 

In the following we present numerical results for the two 
principal molecule orientations, that is, for the molecular 
axis aligned perpendicular (tp — ^) and parallel (ip = 0) 
to the surface. 

First, we consider the capturing of a metal electron 
into the molecular ground state level, which was described 
in Sec. [2] Figure [2] shows the evolution of the ground state 
level's occupation no(t) as calculated from Eq. (29) for the 
case of an AI2O3 surface at a kinetic energy of50mey. 
We clearly see that the process is very efficient, since the 




Fig. 3. Fraction of negative ions versus time for the 
surface-induced decay process (upper panel) and the nat- 
ural decay reaction (lower panel) in front of an AI2O3 sur- 
face at a molecular kinetic energy of 50meV. The solid 
lines represent the parallel molecule orientation, whereas 
the dashed lines denote the perpe ndicu lar orientation. The 
results were calculated from Eq. ( 70b ) . 



molecular vacancy is completely filled within the incom- 
ing branch of the trajectory for both orientations. The effi- 
ciency is higher in parallel orientation. Complete filling of 
the ground state hole is achieved in that case at about t — 
— 1, corresponding to a molecule-surface distance of ap- 
proximately 4.92 ag. In the perpendicular orientation, on 
the other hand, complete filling is only reached at about t = 
which corresponds to the molecule's turning point zq = 
4.42 a B . 

Next, we focus on the decay of a negative ion generated 
by electron capture. To estimate the relative strength of 
the surface-induced and the natural decay channel, respec- 
tively, we treat them separately. Figure [3] shows the time 
evol ution of the fraction of negative ions as calculated from 
Eq. (70b I for both processes in front of an AI2O3 surface 
at £kin = 50 meV. As we see, within the incoming branch 
the portion of negative ions first increases due to the ef- 
ficiency of the electron capture into the ground state and 
the finiteness of the decay rates. At some point within the 
incoming branch, however, the decay process starts to out- 
balance the generating reaction, resulting in a decrease of 
the negative ion fraction. Since the surface-induced d eca y 
process is driven by the molecule-surface interaction pdj ) , 
which vanishes at large distances from the solid, the neg- 
ative ion share saturates to a finite value in this case (see 
upper panel of Fig. [3]). For the parallel orientation the 
saturation level is much lower than for the perpendicu- 
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Fig. 4. Time evolution of the fractions of metastables 
and ground state molecules n g in parallel (solid lines) and 
perpendicular (dashed lines) orientation for an AI2O3 sur- 
face at a molecular kinetic energy of 50 meV. Only the 
natural decay channel was consi dered. Th e curves were 
obtained by employing Eqs. (70a) and (70c), respectively. 




Fig. 5. Spectrum of emitted electrons for an AI2O3 sur- 
face in parallel (solid lines) and perpendicular (dashed 
lines) orientation at t — 00, as calculated from Eq. ( |73[ ). 
The molecule's kinetic energy was fixed to 50meV. Only 
the natural decay process was taken into account. 



lar orientation. This is caused by the antisymmetry of the 
molecular wave functions in the plane normal to the mole- 
cule's axis. In the parallel orientation this antisymmetry is 
broken by the surface, which results in a significantly en- 
larged matrix element, and hence in a higher rate of decay. 
For the natural decay, on the other hand, the decay rate 
is constant, and thus the negative ion fraction eventually 
drops to zero (see lower panel of Fig. |3| . 

An important point to notice is that the portion of neg- 
ative ions is about two orders of magnitude higher in the 
surface-induced decay as compared to the natural decay. 
The natural decay rate must thus be considerably larger 
than the surface-induced decay rate. Hence, the latter can 
be neglected in the present situation. It could, however, 
become important for material combinations that allow 
for projectile turning points closer to the surface. This 
would lead to an increase in the matrix element, which in 
turn results in a higher decay rate. 

To calculate the fractions of metastables n + and ground 
state molecules n g from Eqs. (70a) and (70c I, respectively, 
we focus therefore, in view of the above discussion, ex- 
clusively on the natural decay and neglect the surface- 
induced decay channel. The results are depicted in Fig. [4j 
Note that with the decrease of the fraction of metastables 
the ground state share almost immediately rises by the 
same amount, which is due to the rather large decay rate 
for natural decay. The fractions and n g can be related 
to the occupancies of the molecular levels no and n\ by 
means of Eq. (71 ). 

We now turn to the spectrum of emitted electrons, 



which can be calculated from Eq. ( 73 ) . Figure [5] shows 
the electronic spectrum at t = 00, again for the case of 
AI2O3 and £kin = 50meV, while neglecting the surface- 
induced decay channel. In both molecular orientations the 
spectrum exhibits a strong cut-off for energies below ap- 
proximately 1 eV. This is a direct consequence of the im- 
age potential, trapping low energy electrons close to the 
surface, which was incorporated into our calculation by 



means of the quantity a(s qz ,t) (see Eq. (43)). Following 
the low energy cut-off, the spectra show a strong peak 
at about 1.8 — 1.9 eV, before they slowly fall off towards 
zero for larger energies. Examining the curves we further 
note that the spectrum is larger in the parallel orientation 
than in the perpendicular orientation. This can again be 
attributed to the trapping effect. From Fig. [3] we know 
that in perpendicular geometry the negative ion is gener- 
ated and destroyed closer to the surface than in parallel 
orientation. The image potential and hence the trapping 
effect are thus particularly strong in that case which in 
turn reduces the efficiency of electron emission. 

Finally, we calculate the secondary electron emission 
coefficient 7 e , which is the area beneath curves of the type 
shown in Fig. [5] We specifically consider AI2O3, MgO, 
Si02, and diamond. The electronic parameters for these 
substances are listed in Table [I] We assume the turning 
point formula given in [5] to be valid for all of these ma- 
terials. Hence, we can neglect the surface-induced decay 
process. 

Figure [6] depicts the variation of the secondary elec- 
tron emission coefficient 7 e with the molecule's kinetic en- 
ergy Skin for the aforementioned materials as calculated 



from Eq. ( 73 ) . Leaving MgO aside (for a discussion see be- 
low), the 7 e -cocfficicnts arc on the order of 10 _1 over the 
whole range of kinetic energies. The rather large values 
for 7 e can be attributed to the shape resonance N^ 2 /^) 
which is not only efficiently formed in front of the sur- 
face but also quickly decays thereby releasing an electron. 
Atomic projectiles, for instance, metastable argon, which 
cannot form a metastable negative ion, will lead to much 
smaller j e values. 

The secondary electron emission coefficient we find 
for metastable nitrogen hitting a dielectric surface agrees 
rather well with the 7 e -coefficient which is required for 
a convergent, self-consistent kinetic simulation of DBDs 
taking this particular surface collision process into ac- 
count [3] . In keeping with our initial motivation, the effec- 
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Fig. 6. Secondary electron emission coefficient 7 e for sev- 
eral different dielectric surfaces plotted against the mole- 
cule's kinetic energy Ekin- For each material two curves are 
depicted. The upper corresponds to the parallel molecule 
orientation, whereas the lower identifies the perpendicu- 
lar molecule orientation. All results were obtained from 



Eq. (731. 



tive microscopic model is thus indeed capable to provide 
input data with sufficient accuracy for the kinetic mod- 
elling of gas discharges. 

Let us now analyze Fig. [6] in more detail. The 7 e - 
coefneients increase with decreasing kinetic energy. This 
can be explained by the enlarged molecule-surface interac- 
tion time and the narrowed turning point which lead to a 
more effective filling of the ground state hole at lower pro- 
jectile velocities. Apart from this obvious monotony prop- 
erty, the specific form of the particular graphs is, however, 
not easily understood. The reason for this is the complex 
dependence of the secondary electron emission coefficient 
on the surface's band structure and its static dielectric 
constant. For the materials we considered a few general 
remarks, however, can be given. 

First, we note that the emission coefficient is sensitive 
to the efficiency of the electron capture into the molecular 
hole level. A high efficiency results in negative ions getting 
produced at large distances from the surface. The subse- 
quent electron emission will thus also take place at large 
distances, because of the large natural decay rate of the 
ions. This is beneficial to the emission coefficient, as the 
image potential trapping is less severe at larger distances. 
The efficiency of the electron capture in turn is determined 
by the alignment of the solid's valence band to the molec- 
ular hole level. Since the filling of the molecular vacancy 
is a resonant tunneling process, the vacancy level Eq has 
to be degenerate to occupied states within the surface in 
order to allow for an efficient generation of negative ions. 
In addition, the electron capture is more productive if the 
degeneracy appears at higher energies within the valence 
band. This is due to the fact that at higher energies the 
wave functions of the electrons within the solid have a 
larger extension outside the solid. Thus their overlap with 
the molecular wave function is increased, which directly 
influences the matrix element and thus the efficiency of 
the electron capture. For MgO both conditions, degener- 



acy and higher energies within the valence band, are vi- 
olated since the molecular hole level Sq shifts downwards 
out of the valence band at small distances (see Fig. [IJ . As 
a result MgO has the smallest emission coefficient of the 
materials under consideration. 

Secondly, a higher dielectric constant seems to have an 
adverse effect on the value of j e , since for the most part the 
emission coefficients for Si02 and diamond, both having 
lower dielectric constants, are higher than the emission 
coefficient of AI2O3 which possesses a higher dielectric 
constant (see Table [I]). This effect can probably be at- 
tributed to the amount of image potential trapping of slow 
electrons, which is proportional to [e b r — l)/(£r + 1) (see 
Eqs. ^ and (43)). This factor increases with e h r . Thus, 
the image potential and with it the effect of trapping is 
enlarged for higher dielectric constants. 



5 Conclusion 

In this work we described a generic effective microscopic 
model for secondary electron emission due to collision of 
metastable molecules with dielectric surfaces. As an illus- 
tration we applied the model to N2( 3 I?+) hitting a di- 
electric surface, focusing on the RCT channel |2]) which 
is particularly strong because this molecule forms a shape 
resonance upon electron capture. In addition, the com- 
peting Auger process ([!]) is energetically blocked for the 
dielectrics commonly used in gas discharge physics and 
needs thus not to be considered. 

Our model depends only on a small number of mate- 
rial parameters. As far as the surface is concerned, these 
parameters include the electron affinity e a of the solid, its 
band gap e g , and the width of its valence band Aey as well 
as its static dielectric constant e\. Any surface that can 
be parameterized by such a set can in general be treated 
within our model. With a few modifications the model can 
also be applied to metallic surfaces and other molecular 
and atomic projectiles. 

We presented numerical results for AI2O3, MgO, Si02, 
and diamond surfaces. The electronic parameters of the 
surface in question are crucial for the overall efficiency 
of the RCT process. The generation of the shape reso- 
nance N^^TTg) upon electron capture into the lower ion- 
ization level of N^~( 2 77 g ) is efficient only when the molec- 
ular level stays well inside the valence band. Concerning 
the decay of the temporary negative ion we found the 
surface-induced process to be about two orders of mag- 
nitude weaker than the auto-decay. The surface-induced 
decay channel can thus be neglected for the situations 
we investigated. We found the spectrum of the emitted 
electron to be strongly peaked at about 1.8 — 1.9 eV. In 
a gas discharge, however, the auto-decay occurs close to 
the plasma wall where the sheath potential is large. The 
emitted electron will thus immediately gain a kinetic en- 
ergy which is substantially larger than the energy resulting 
from auto-decay. 

Calculations of the secondary electron emission co- 
efficient 7 e for different surface materials finally yielded 
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7 e -values from 0.02 to 0.2 for collision energies ranging 
from 0.1 fieV to 1 eV. Particularly for thermal energies 
the 7 e -values we obtain are of the order of 10 _1 and thus 
coincide with the values deduced from kinetic simulations 
of dielectric barrier discharges which take this particular 
secondary electron emission process into account. Hence, 
despite its simplicity effective microscopic modelling can 
capture the essential physics of secondary electron emis- 
sion from surfaces due to impacting mctastablc molecules 
and is thus capable to fill some of the gaps in the molecule- 
surface collisions' reference data which are required for the 
kinetic modeling of gas discharges. 
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